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ABSTRACT 

Using data from the SCUBA Legacy Catalogue (850 ;^m) and Spitzer Space Telescope 
(3.6 - 70 ^m), we explore dense cores in the Ophiuchus, Taurus, Perseus, Serpens, and 
Orion molecular clouds. We develop a new method to discriminate submillimeter cores 
found by SCUBA as starless or protostellar, using point source photometry from Spitzer 
wide field surveys. First, we identify infrared sources with red colors associated with 
embedded young stellar objects (YSOs). Second, we compare the positions of these 
YSO-candidates to our submillimeter cores. With these identifications, we construct 
new, self-consistent starless and protostellar core mass functions (CMFs) for the five 
clouds. We find best fit slopes to the high-mass end of the CMFs of —1.26 it 0.20, 
-1.22 ± 0.06, -0.95 ± 0.20, and -1.67 ± 0.72 for Ophiuchus, Taurus, Perseus, and 
Orion, respectively. Broadly, these slopes are each consistent with the —1.35 power-law 
slope of the Salpeter initial mass function (IMF) at higher masses, but suggest some 
differences. We examine a variety of trends between these CMF shapes and their parent 
cloud properties, potentially finding a correlation between the high-mass slope and core 
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temperature. We also find a trend between core mass and effective size, but we are very 
limited by sensitivity. We make similar comparisons between core mass and size with 
visual extinction (for Ay > 3) and find no obvious trends. We also predict the numbers 
and mass distributions of cores that future surveys with SCUBA-2 may detect in each 
of these clouds. 



1. Introduction 

The origin of stellar mass is not well understood. We know stars form in very cold and dense 
regions of gas and dust, deeply embedded within molecular clouds, but the dominant processes 
that guide the formation of small scale structure and the collapse of those structures into stars 
is unclear. Wc define zones of relatively high density on the small scales in molecular clouds as 
"cores," which have masses ~ Mq within a radius of 0.07 pc and would form a single star or a 
stellar system with a few stars (Di Francesco et al. 2007). Thus, molecular cloud cores represent the 
cumulative result of all physical processes that organize cloud mass into relatively high densities 
on small scales. 

The masses of molecular cloud cores are important probes of the initial conditions of star 
formation, and the relationships between these cores and any stellar products may be key to 
understanding the origin of stellar mass. Previous studies of these populations (for examples, 
see Motte et al. 1998; Nutter & Ward-Thompson 2007; Ward-Thompson ct al. 2007a; Enoch 
et al. 2008) have revealed that core mass functions (CMFs) roughly follow lognormal shapes that 
closely resemble the observed stellar initial mass function (IMF). In particular, these studies have 
shown great similarity between the high-mass slopes of CMFs and the high mass slope of the IMF. 
The CMF, however, may still depend on physical conditions within the parent cloud. Indeed, 
the relationship between the CMF and IMF is likely complex and a variety of factors, such as 
fragmentation (e.g., Dobbs ct al. 2005), competitive accretion (e.g., Bonnell et al. 2004), turbulence 
(e.g., Elmegreen 2002), magnetic fields (e.g., Shu et al. 2004), and radiative feedback (e.g., Offner 
et al. 2009) may be relevant. 

The continuity between CMFs and the stellar IMF suggests that stellar mass is related to 
how material in molecular clouds is first collected into stellar precursors. Unfortunately, previous 
studies generally examined small regions within larger clouds (e.g., the L1688 region of Ophiuchus 
in Motte et al. 1998) and over different wavelengths (e.g., 850 fim, 1.1 mm). As well, each group 
had developed a different set of conditions to identify small scale structure in molecular clouds, 
thereby these studies follow different numbers of cores, even for the same region (e.g., for Perseus, 
Hatchell et al. 2005; Kirk et al. 2006). In spite of the inconsistent analyses, these previous studies 
seem to produce CMFs that resemble the IMF shape. Nevertheless, it is difficult to compare CMFs 
from different studies and test for differences between clouds. It would be advantageous to look at 
core populations in various clouds with different physical conditions using a self-consistent analysis. 
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In addition, observations of dense cores over the last decade have revealed populations of 
cores with and without embedded young stars (Di Francesco et al. 2007). Cores that contain 
a young stellar object (YSO) will have lost some of the surrounding material to outflows or to 
accretion onto the central body (Myers 2008). Also, their intrinsic temperatures may differ, biasing 
estimates of their masses. Thus, to obtain an accurate CMF, starless cores must be differentiated 
from protostellar cores, ie., those containing YSOs. 

Distinguishing between starless and protostellar cores depends on detecting an infrared source 
embedded in the dense material. Deeply embedded sources can be very faint, however, making 
such detections difHcult. As well, many observed infrared sources are background galaxies or bright 
giant stars in our Galaxy, which must be identified as external to the cloud. Therefore, multi-band 
mid- to far-infrared photometry has emerged as an effective means for identifying and classifying 
embedded YSOs from external sources of infrared emission (e.g., Allen et al. 2004; Megeath et al. 
2004; Harvey et al. 2006; Gutermuth et al. 2008; Megeath et al. 2009; Evans et al. 2009; RebuU 
et al. 2009). 

Already, numerous studies have developed techniques to identify protostellar core populations 
from starless core populations in molecular clouds (e.g., J0rgensen et al. 2007; Enoch et al. 2009). 
These techniques, however, employ very different constraints in infrared colors, magnitudes, or 
locations of YSO candidates with respect to the cores. As such, the CMFs produced by each 
technique will contain different biases, making comparisons between these studies unclear. Ideally, 
we would want to use the same constraints to compare CMFs over a large sample of molecular 
clouds to explore differences and similarities in core populations with environment. 

Here, we use 850 jim data from the Submillimeter Common-User Bolometer Array (SCUBA) 
and just-released Spitzer data with a new self-consistent method of classifying cores to identify 
starless and protostellar cores in five clouds (Ophiuchus, Taurus, Perseus, Serpens, and Orion). 
With this discrimination, we produce consistent CMFs across the five clouds and thus, we can 
directly compare these distributions. In §2, we discuss our sample choices, including the target 
clouds, the infrared and submillimeter data used in this study, and the individual core populations 
and selection criteria. In §3, we discuss our new core classification technique and compare our 
results with two other previously developed methods. In §4, we examine the CMFs produced from 
our own classification method, and we compare these to standard formulations of the IMF. Also, we 
examine trends in the CMFs between the clouds, compare core properties with their surrounding 
environments, and make predictions as to what forthcoming instruments will detect. In §5, we 
summarize our results. 
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2. Clouds 

2.1. Cloud Properties 

Our analysis focused on the Ophiuchus, Taurus, Perseus, Serpens, and Orion molecular clouds. 
These clouds were chosen due to a wealth of available data (e.g., millimeter, infrared) and their 
close proximity (all < 500 pc), resulting in maps with good linear resolution. Such small scale 
observations are necessary to resolve cores as well as to probe the physical properties and structure 
inside cores (Ward-Thompson et al. 2007a). The five clouds studied here represent a variety of 
physical environments. For example, the Taurus cloud is undergoing only low-mass star formation 
(Hartmann 2000) whereas the Orion region contains active regions of massive star formation and 
is found within a large OB association (Peterson & Megeath 2008). We found total cloud masses 
using our visual extinction maps (see §2.3.3 and Figure 1). Assuming that 1 magnitude of extinction 
corresponds to 10^^ molecules of hydrogen gas (1 Ay = 10^^ H2 cm~^) and a mean molecular weight 
of 2.33, we measured the mass of the cloud via, 

M(M,) = 56.4 (^)'(^)'i:^v,< (1) 

where 6 is the angular size of the extinction map pixels, d is the cloud distance (listed in Table 
1), and Av,i is the extinction in a given pixel element. Table 1 gives the mass measurements we 
obtained for each cloud and the area over which we measured the mass. We constrained our cloud 
mass and size measurements to the regions with Ay > 1. For Serpens, the cloud is located close 
to the Galactic Plane, making its boundary less certain. Thus, we are likely including extinction 
from the nearby Auriga cloud and over-estimating its mass. 

2.2. Core Properties 

For each cloud, we assumed the core dust temperatures, T^i, and the 850 fim opacities, kssO) 
were constant. Temperatures and opacities of core dust probably deviate within a given region due 
to different circumstances, such as proximity to an embedded cluster. Table 2 lists the assumed 
values for and Table 1 gives the assumed values of distance for each cloud. 

We note that the temperatures listed in Table 2 were obtained via different techniques. Friesen 
et al. (2009) and Rosolowsky et al. (2008) derived kinetic temperatures, Tk, of dense gas in Ophi- 
uchus and Perseus, respectively, using ammonia hyperfine structure lines, and we assume these 
temperatures are similar to since the densities in the cores are expected to be high (Goldsmith 
2001). For the L1544 region in Taurus, Andre et al. (2000) used spectral energy distribution (SED) 
fitting from ISO, SCUBA, and IRAM observations to obtain = 13 K. As part of the COM- 
PLETE survey, Schnee et al. (2005) used the 60//m/100/Lim flux density ratio for Serpens (see their 
Figure 5) to estimate its dust temperature. For Orion, Wilson et al. (1999) observed ammonia 
and derived rotational temperatures of ~ 20 K. As with Ophiuchus and Perseus, we assume the 
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rotational temperature is similar to the dust temperature since the gas traced by ammonia is dense 
(~ 10'' cm-3). 

Clearly a common origin of would be preferable, but note that the masses of all cores in the 
CMFs will scale with and to first order the CMF shape will not depend on (see §4 for further 
discussion; we assume a 30% error in to derive uncertainties in our masses and CMFs). Self- 
consistent determinations of for cores in these clouds will soon be possible through SED fitting of 
75 — 500 /xm data from the Herschel Gould Belt Survey (Andre k. Saraceno 2005; Ward-Thompson 
et al. 2007b). 

2.3. Data 

We obtained our data from large-scale surveys and included wavelengths from the submillimeter 
(850 jJLva) to the infrared (~ 2 jim to 70 jim). We used the submillimeter data to probe the densest 
regions of each cloud and the infrared data to identify embedded protostars through emission and 
characterize the extended cloud structure through extinction. We discuss each of these data sets 
in turn below. 

2.3.1. SCUBA Maps 

We obtained our submillimeter data from the SCUBA Legacy Catalogue (SLC; see Di Francesco 
et al. 2008 for more details), which is an archive of all 850 /xm and 450 ^m. mapping observations 
from SCUBA on the James Clerk Maxwell Telescope (JCMT). These observations generally focused 
on the high extinction regions within each cloud, thus the SLC data are patchy. The SLC includes 
two sub-catalogues: the Fundamental Catalogue, which contains only objects identified from data 
with high quality atmospheric corrections (consisting of ~ 78 % of map data with an arcal coverage 
of ~ 19.6 deg^), and the Extended Catalogue, which includes all the data regardless of quality (areal 
coverage of ~ 29.3 deg^). For consistency, the SLC used the algorithms developed by Johnstone 
et al. (2000a) to similarly reduce these data. As well, the SLC used the 2D Clumpfind algorithm 
(Williams et al. 1994) to identify structures in the continuum emission. Clumpfind identifies flux 
peaks over a certain noise level and then uses flux contours (down to a minimum threshold level) 
to assign boundaries to the cores and measure total fluxes (see Di Francesco et al. 2008). For a 
discussion on the biases associated with different clump-finding algorithms, see Curtis k. Richer 
(2009). 

To identify our submillimeter cores, we only used the 850 //m data, since the 450 //m obser- 
vations have a greater absolute flux uncertainty than the 850 ^im. data by over a factor of two (Di 
Francesco et al. 2008) . Since we are more interested in accurate core fluxes than wide areal coverage 
of the clouds, we drew our samples from the Fundamental Catalogue. Table 3 lists the cloud areal 
coverage observed by each of the surveys, where the second column indicates the area mapped in 
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the Fundamental Catalogue. 

Initially, the SLC set a Clumpfind threshold level equal to a factor 3 above the local noise level 
in a given map. For each of these objects, a second set of flux and size measurements, known as the 
"alternative" flux and radius, were calculated assuming a common threshold of 90 mJy beam~^, 
which is a factor 3 larger than the typical 850 fim noise level of all SLC maps, 30 mJy bcam^^ 
(Di Francesco et al. 2008). For our core properties, we considered the "alternative" flux and radius 
measurements given in the SLC to provide a consistent mass sensitivity for all the cores in a given 
cloud. Generally, the two flux and size measurements were quite similar. 

2.3.2. Spitzer Space Telescope Maps 

The Spitzer Space Telescope (Werner et al. 2004) has two instruments that observe mid- and 
far-infrared wavelengths: IRAC (Infrared Array Camera; Fazio et al. 2004) at 3.6, 4.5, 5.8, and 
8.0 /im and MIPS (Multiband Imaging Photometer for Spitzer; Rieke et al. 2004) at 24, 70, and 
160 /um. With the high infrared sensitivity provided by Spitzer, these cameras provided excellent 
data for determining the presence of a protostar within highly extincted regions like cores. We 
used "Molecular Cores to Planet Forming Disks" (c2d; Evans et al. 2003, 2009) Legacy Project 
data for the Ophiuchus, Perseus, and Serpens molecular clouds (see J0rgensen et al. 2006; Harvey 
et al. 2006; Rebull et al. 2007; Harvey et al. 2007a,b; Padgett et al. 2008). The Taurus data came 
from Guest Observer (GO) observations for the Taurus Legacy Project (Padgett et al. 2010 in 
preparation), and the Orion data arc from Guaranteed Time Observations (GTO) for the Spitzer 
Orion Survey (Megeath et al. 2010 in preparation). 

For the c2d clouds, Taurus, and most Orion observations, IRAC and MIPS observations were 
conducted over two epochs to identify and remove asteroids. Because the observations were con- 
ducted over two epochs at least 6 hours apart, even distant asteroids can be removed from the 
images and subsequent catalogues. Asteroids are prominent at 24 /xm and 8 /xm and recognizable 
in the other bands. Taurus, in particular, is located towards the ecliptic plane, and about half of 
the point sources detect at 24 /xm in Taurus are asteroids (see Padgett et al. 2010 in preparation). 

For the c2d catalogue, the MIPS integration times were 3 seconds per sky pointing, with a 
given position observed 5 times for a total of 15 seconds. While the 24 ^m and 70 /im bands covered 
roughly the same area of the sky, the 24 ^m scans had longer total integration times (30 seconds) 
than the 70 //m scans (15 seconds; Young et al. 2005). Observations at 160 //m were not included in 
the final c2d catalogues, since these data were affected by saturation and a large beam size (Evans 
et al. 2009). Similar to MIPS, IRAC observed each cloud twice, with the two epochs separated by 
about 6.5 hours. The integration time per pointing with IRAC was 12 seconds (Porras et al. 2007). 
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The Taurus observations were conducted and described by Padgett et al. (2008^, 2010 in 
preparation), and consisted of IRAC and MIPS observations over a wide area (see Table 3). While 
the MIPS observations were obtained in the same fashion as the c2d maps, the IRAC maps sacrificed 
redundancy for area and consisted of only two (12-second) frames per position. At the distance of 
Taurus (140 pc), however, these shallow IRAC maps should still detect legitimate YSOs down to 
very low masses, unless these sources are exceptionally embedded. 

For Orion, the MIPS observations are a combination of several different surveys with different 
observing modes. The MIPS observations consisted of both fast (17" s~^) and medium (6.5" s~^) 
scanning modes with total integration times ranging from 30 - 40 seconds. Similar to the c2d IRAC 
observations, several regions in Orion were observed with IRAC in high dynamic range mode, 
though some were observed only once in a faster frame time. Only the IRAC observations for 
Orion were made over two epochs to minimize artifacts and maximize four band coverage. For 
more information, see Megeath et al. (2010, in preparation). 

The areal coverage of our five clouds by Spitzer are listed in column 3 of Table 3. Generally, 
MIPS and IRAC scans did not cover identical regions, though there was significant overlap. Due to 
these different areal coverages and different effective sensitivities among the instruments, there are 
a number of sources that were detected at only a few wavelengths. For objects that were identified 
in the IRAC or MIPS maps but without valid photometry at a particular band, the c2d team 
"band-filled" the missing band by using a wavelength-appropriate point spread function (PSF) 
located at the expected position of the object. In cases where the pixels are compromised (e.g., 
from saturation or cosmic ray contamination) or the object is undetected, however, this process 
can result in negative band-filled fluxes. For more information, see the c2d Delivery Document, 
Evans et al. (2007). The Taurus and Orion data are not band-filled. 

2.3.3. 2MASS Extinction Maps 

Infrared observations can probe the total line-of-sight column density of dust by tracing the 

color of stars located beyond a molecular cloud that are reddened due to high dust extinction (e.g., 
Lada et al. 1994). Extinction maps for each of our five clouds were created using archived 2-Micron 
All Sky Survey (2MASS) catalogues of point sources in a similar manner described in Lada et al. 
(1994), Lombardi k Alves (2001), Cambresy et al. (2002), and Schneider et al. (2006). For more 
details, see Schneider et al. (2006). 

First, individual extinction values were obtained from a weighted average of the J-H and H- 
K colors of individual stars, assuming the average intrinsic colors were (J-H)o = 0.45 it 0.15 and 
(H-K)o = 0.12 lb 0.05 as derived from stellar population models and typical dispersions of galactic 
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stars measured using the Besancon Galactic modeP by Robin et al. (2003). Second, the Besancon 
stehar models were used to derive average densities of foreground stars in the 2MASS bands at the 
distance of each cloud. For clouds near the Galactic plane or at large distances, the foreground star 
populations can be significant. The expected numbers of foreground stars were removed using the 
bluest sources (see Cambresy et al. 2002 for more details). Finally, visual extinction was derived 
from averaging the individual Ay measurements for stars within a Gaussian beam. This averaging 
excluded the bluest sources and was adapted so that > 10 stars contribute to the extinction element. 

2.3.4- Regions with Complete Submillimeter - Infrared Coverage 

As an all-sky survey, 2MASS data naturally encompassed the full extent of our 5 clouds. The 
Spitzer coverages of these clouds were quite large, but generally restricted to areas of Ay > 3. 
SCUBA was used to map wide regions only rarely (see Johnstone et al. 2004, Hatchcll et al. 2005, 
Kirk et al. 2006) given its limited sensitivity and the observations typically focused on regions of 
known star formation. 

Figure 1 compares the areal coverage from 2MASS (background) with the areal coverage from 
the Fundamental SLC (contours). Also included are boundaries to distinguish the Scorpius and 
Monoceros clouds from the Ophiuchus and Orion clouds in their respective maps. The boundaries 
for Serpens are not clear and much of our extinction map is associated with the Auriga system. 
From Figure 1, we see that much of these clouds remain unmapped in the submillimeter, including 
some regions of high extinction. Table 3 lists the cloud area observed by each of the surveys. 

2.4. Identification of Cores Within Clouds 

Using the Fundamental SLC we generated lists of core candidates for each of our clouds. Since 
our extinction maps extend over entire clouds, we used these data to determine where the core 
candidates were located within the large-scale structure of their parent clouds. We did this by 
identifying the extinction level at the position of each of our core candidates. For this study, we 
required that cores be located in a cloud region of Ay > 3 to ensure the cores arc dense. Generally, 
SCUBA observed higher extinction regions, so this constraint did not greatly affect the number of 
core candidates in our sample. 

Also, we visually inspected our core candidates to remove those that were likely artifacts 
of flat-fielding or which appeared too diffuse to be dense cores. Finally, to ensure we had good 
detections, we removed all objects with 850 //m peak fluxes less than 5a (5 x 30 mJy beam~^ = 150 
mJy beam"-*^), where a is the average noise level of 30 mJy beam~^. Table 4 summarizes all the 
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cuts made to the core candidates extracted from the SLC Fundamental Catalogue. For Serpens, 
only one relatively small region of ~ 1 pc^ was observed with SCUBA, resulting in far fewer core 
numbers relative to the other clouds. 

3. Separating Starless and Protostellar Cores: A New Classification Technique 

3.1. Summary of Previous Approaches 

In compiling samples of cores observed in each cloud, we did not distinguish between differ- 
ent stages in core evolution, such as those that are starless or those that are protostellar. These 
populations must be separated to explore properly the relationship between the IMF and CMF 
(Ward-Thompson et al. 2007b). Previous efforts have attempted to separate the starless and pro- 
tostellar core populations in molecular clouds by using infrared observations with (sub)millimeter 
continuum data. Slightly different methodologies were used, however, to accomplish this separation. 

There are many possible approaches to identifying protostellar cores. For example, both the 
J0rgensen method (see J0rgensen et al. 2006, 2007, 2008; hereafter J06, J07 and JOS, respectively) 
and the Enoch method (see Enoch et al. 2008, 2009; hereafter EOS and E09, respectively) considered 
a (sub) millimeter core to be protostellar when a YSO-like infrared source was in close proximity to 
the dense core. These two methods, however, differed in their identification of YSO-like infrared 
sources and in their definition of "close proximity". We discuss differences between these two 
techniques in §4.1.3. Aside from these methods, other approaches included using the SEDs of point 
sources to identify infrared excesses, calculate spectral indices or measure bolometric temperatures 
(e.g., Hatchell et al. 2007; Winston et al. 2007; Evans et al. 2009) and using infrared colors to 
identify and classify YSOs (e.g., Allen et al. 2004; Megeath et al. 2004; Harvey et al. 2006; Kirk 
et al. 2009; Gutermuth et al. 2008). As a hybrid approach, Megeath et al. (2009) used colors to 
measure SED slopes and identify protostars. 

Which approach should be used? All of these methods are subject to uncertainties from un- 
known reddening levels and possible chance coincidences. Overall, SEDs can accurately assess 
whether cores are starless or protostellar (ie., through the bolometric temperature or indicators 
such as infrared excesses from disks), but this requires a wealth of high-resolution data at a variety 
of wavelengths, which is observationally expensive. In addition, resolution can vary over differ- 
ent wavelengths, thus complicating how objects detected at different wavelengths are associated. 
Combining color and spatial co-location criteria would be least biased to particular data sets such 
as those with high resolution or large spectral coverage. Accordingly, we describe our new core 
identification system in the next sections. 
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3.2. New Color Criteria 

To produce a starless CMF, one must identify and omit cores that have lost some of their 
surrounding envelope from outflows or from accretion onto a protostar. This task is complicated 
because common sources of infrared emission include galaxies and background stars along the line 
of sight. Fortunately, criteria have been developed for distinguishing YSOs from star forming 
galaxies, active galactic nuclei (AGNs), and giant stars (e.g., Harvey ct al. 2006; E09; Gutermuth 
et al. 2008; Megeath et al. 2009; Kirk et al. 2009; Rebull et al. 2009). In general, extragalactic 
sources are often very faint or have unique colors, such as particular spectral signatures due to 
strong polycyclic aromatic hydrocarbon (PAH) emission associated with HH regions, causing an 
increased infrared excess at 5.8 fim and 8.0 /xm. PAH emission features are not readily detected 
towards embedded protostars, ensuring that extragalactic sources are identified and not protostellar 
sources (van Dishoeck & van der Tak 2000; Geers et al. 2009). Furthermore, Spitzer observations 
can distinguish between objects with infrared excesses from heated dust grains and from pure 
photospheric emission of reddened background stars. Red stellar sources, such as AGB stars, are 
likely to have flat spectra in the infrared, thereby distinguishing them from protostellar sources, 
which peak in the mid- to far-infrared (Harvey et al. 2006). 

Young protostars still embedded in dense cores should peak in the far-infrared due to the 
absorption and reprocessing of light to longer wavelengths by envelopes (a characteristic envelope 
temperature is T ~ 20 K). Several recent studies have identified these sources using IRAC and 
MIPS colors (e.g., Harvey et al. 2006; J07; Megeath ct al. 2009) or 2MASS colors (e.g., HatcheU 
et al. 2007; Kirk et al. 2009). Given the scatter and overlap of various objects in color or magnitude, 
however, there will be some objects that are not selected and conversely, not all contaminants will 
be removed. In addition, excesses in some bands could be the result of different physical processes. 
For example, protostellar outflows can interact with the surrounding molecular cloud and result in 
strong emission in the IRAC bands due to shocked H2 gas (Gutermuth et al. 2008). 

To identify embedded YSOs from the millions of compact infrared sources detected by Spitzer, 
we chose specific infrared sources based on the quality of the detection and the source colors. Since 
embedded protostars peak in the far-infrared, we required that our YSOs have strong detections at 
24 ^m and 70 fim and rising-red colors. The red colors will exclude stellar sources, which have flat 
colors in the infrared regime, but not extragalactic sources. Extragalactic contamination must be 
excluded separately. Our complete color criteria (CC) for identifying protostellar objects is listed 
below: 

CCl. source 24 fim or 70 /xm flux has a S/N > 5 or if the 70 /xm data was poor, a S/N > 5 in 
all four IRAC bands, and 

CC2. either 24 /tm or 70 /tm detections are real (not band-filled limits), and 

CCS. source colors are dissimilar to those of star-forming galaxies (see Gutermuth et al. 2008), 
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CC4a. if detected at 24 //m, [8.0] - [24] > 2.25 and [3.6] - [5.8] > -0.28([8.0] - [24]) + 1.88, or 
CC4b. if not detected in 24 //m, [3.6] - [5.8] > 1.25 and [4.5] - [8.0] > 1.4 

Since embedded protostars should emit strongly in the far-infrared, CCl selects infrared sources 
that were well detected at 24 /xm or 70 fim. MIPS, however, has a lower sensitivity than IRAC, so 
we may miss sources relying on strong MIPS emission. Indeed, for Orion and Taurus, there were 
relatively few unique 70 ;um detections with respect to unique 24 nm detections. For these clouds, 
we modified CCl to identify sources based on S/N > 5 either at 24 nm or at all four IRAC bands, 
rather than using 70 /xm. For CC2, we followed E09 and excluded sources that were band-filled (see 
§2.3.2) at either 24 /itm or 70 /xm. Since MIPS observations have a lower resolution than the IRAC 
bands, bandfilling at 24 //m or 70 /im could result in contamination from the wings of nearby bright 
sources. Such observations make the MIPS fluxes uncertain. CC3 excludes star-forming galaxies 
based on the prescription developed by Gutermuth et al. (2008), which selects for a growing infrared 
excess at 5.8 //m and 8.0 /xm due to strong PAH emission. 

The colors given in CC4a and CC4b were determined using the results from the c2d survey 
(Evans et al. 2009). Evans et al. (2009) produced a very large sample of YSOs (1024 over 5 clouds) 
classified into Class O/I, Flat, Class II, and Class III using the source SEDs. In their Figure 11, 
they plotted their sample in [3.6] - [5.8] and [8.0] - [24] color-color space (with different symbols 
for the different classes). Since the different classes clustered together into specific regions of the 
color-color plot, we can use the [3.6] - [5.8] and [8.0] - [24] infrared colors to identify the degree to 
which our infrared sources are embedded. We wish to identify all young protostars still embedded 
in a dusty envelope, which includes Class 0, Class I, and Flat spectral source types (J08; Myers 
2008). As such, we identified the [3.6] - [5.8] and [8.0] - [24] colors from Figure 11 of Evans et al. 
(2009) to include these spectral classes (CC4a). For objects without 24 /xm detections (ie., due to 
the lower sensitivity), we used CC4b which adopts the [4.5] - [8.0] color criteria from Harvey et al. 
(2007a) and the [3.6] - [5.8] color when [8.0] - [24] = 2.25 (from CC4a). 

We also attempted to remove broad-line AGN contaminants using colors outlined in Gutermuth 
et al. (2008) but found that known young protostars were frequently labeled as AGN by these 
criteria. For example, using this prescription and the Perseus infrared sources, the Class objects 
HH 211, IC 348 MM, and N1333 IRAS 4B (Froebrich 2005) could be misidentified as AGN. Several 
authors (e.g., Winston et al. 2007; Gutermuth et al. 2008; Megeath et al. 2009) apply magnitude 
cuts to reduce the contamination from extragalactic sources. These cuts can also remove faint 
emission from deeply embedded protostars. Extragalactic sources are randomly distributed in the 
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sky. By imposing a condition that our YSO-candidates are coincident with cores (§3.3), we expect 
a low probability of extragalactic contamination, and we do not apply a magnitude limit in our 
selection of YSO-candidates. Studies in J07 and JOS suggest that the likelihood of coincidences 
within 15" of a core by a random distribution of galaxies is only ~ 1%. For the most part, we will 
be considering larger areas than 15" (see §3.3), such that the chance of random coincidences for 
our data could be ~ 4%, though this likelihood is probably lower since we removed sources that 
are likely to be be considered star-forming galaxies using CCS. 

Table 5 lists the number of infrared sources that met each of our color criteria. Our final 
source list may still include some galaxy candidates and other suspect objects, which are unlikely 
to be coincident with a SCUBA core. Figure 2 illustrates the colors for all objects that passed 
CCl and CC2. The solid lines illustrate the CC4a limits. The small light grey circles represent the 
infrared sources found at locations with Sgso < 90 mJy beam^^ that were cut from our sample. The 
crosses represent the infrared sources identified as star-forming galaxies with PAH emission (CC3). 
The diamonds represent our final set of infrared sources. Prom Figure 2, it is clear that some 
star-forming galaxies have [3.6] — [5.8] and [8.0] — [24] colors that match our criteria for protostellar 
cores. This stresses the importance of removing extragalactic sources as well as the difficulty in 
identifying them. 

Another complication is excluding more evolved YSOs (such as Class II soTirccs) which are 
not randomly distributed but found within molecular clouds. Megeath et al. (2009) calculated 
colors using power-law SEDs which had slopes that reflected the different classes of YSOs. For 
example, they find [8.0] — [24] > 3.201 corresponds to Plat (—0.3 < a < 0.3) spectral sources and 
[8.0] — [24] > 3.917 corresponds to rising {a > 0.3) spectral sources (ie.. Class or Class I). In 
comparison, we use a bluer value of [8.0] — [24] > 2.25, which suggests that we may be including 
some Class II sources in our list. We suspect, however, that this effect is minimal since most of 
our sources have [8.0] — [24] > 3.201 (see Figure 2) and most of the sources with bluer [8.0] — [24] 
colors are very red in [3.6] — [5.8]. 

3.3. New Coincidence Criteria 

One of the main conclusions of JOT was that protostars are found very close to the peaks of 
submillimeter emission (see their Figure 2) , indicating that the motion of protostars relative to their 
natal cores is likely small. Therefore, it may be superior to consider a given core to be protostellar 
only if a YSO candidate is found relatively close to the peak intensity. It is certainly reasonable 

to assume, for the most part, that protostars will form in the highly concentrated parts of cores 
which are generally associated with the peak submillimeter flux. This could become problematic 
in instances when many small cores are blended together, at which point the peak value could be 
off-centre from all of the individual cores. 

One important consideration is that cores are typically not circular in projection, but can be 
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very irregular in shape. Thus, a circular approximation of the core extent by using its "effective" 
radius, which can be quite large (~ 50"), or a fixed angular size (ie., 15" as used in the J0rgensen 
method) could probe beyond the respective boundaries of the cores, such as when a core is very 
elongated. In using the effective radius, a larger core has a greater intrinsic chance of coincidence 
with a nearby infrared source. Fixed angular sizes, however, cover different physical scales at 
different distances. 

To ensure that the observed size and shape of the core is considered, we suggest a scheme where 
the object location is compared to a percentage of the difference between the peak 850 fxm intensity 
and the boundary intensity (90 mJy beam~^). Figure 3 demonstrates such an implementation, 
where the search area for coincidence is given by the shaded region. We parameterized this shaded 
region using Equation 2, namely; 

Si limit ~ ^peak ^{.Speak Abound) j (2) 

where Sumit is the minimum associated 850 /xm flux for a YSO candidate to be considered associated 
with a SCUBA core. Speak is the peak 850 fim intensity of the core, Sbound is the threshold intensity 
of the core boundary (90 mJy beam"-*" for the alternative flux), and e is a parameter for scaling 
Siimit between Sp^ak ^ind Shound- Thus, e parameterizes the proximity of a YSO-selectcd infrared 
source to the peak flux of a submillimeter core. With this coincidence criteria, we consider both 
the core size and the flux distribution (core shape). 

We chose a value of e = 0.75 after testing several different limits. Taking e = 0.75, while 
technically arbitrary, corresponds to a fractional size where the number of positive coincidences 
does not rise rapidly and this value is also distinct enough from the boundary (when e = 1) to limit 
the effect from uncertainties in the extent of the core. As such, we considered only cores with at 
least one unique YSO candidate interior to the contour for e = 0.75 to be protostellar. Identifying 
unique YSO candidates within the contour defined by e = 0.75 was not always clear, however. For 
example, some protostellar core candidates were located in filamentary structures or in crowded 
regions which made the surrounding SCUBA flux levels much higher. As such, there were several 
cases where a YSO candidate fell within the e = 0.75 contour of two cores. For YSO candidates 
associated with multiple cores, we assigned the infrared-selected source to the nearest core. 



3.4. Method Comparison 

Table 6 compares the numbers of protostellar cores identified using our criteria (see §3.2 and 
§3.3), with results from applying the J0rgensen and Enoch core classification methods to our SLC- 
derived core lists. The J0rgensen method identified cores as protostellar if a MIPS source was found 
within 15" (1 FWHM of the unsmoothed, non-SLC SCUBA beam) of the core centre, regardless of 
position angle. They also considered infrared sources of suitably red colors and cores of high central 
concentration, which we do not include for our implementation of the J0rgensen method since the 
former is not constrained to submillimeter cores (our aim is to identify starless and protostellar 
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cores) and the latter was shown to be a poor indicator for core evolution (see JOS). Conversely, the 
Enoch method focused on infrared sources with certain properties (ie., certain magnitude levels, 
c2d class designations, and spectral indices) that were located within an intrinsic core FWHM (a 
geometrical mean of the core's major and minor axes), of the core centre, regardless of position 
angle. Since Taurus and Orion were not part of the c2d survey, we could not use the Enoch method, 
which selected sources with favorable c2d class designations^, however, we were able to apply the 
J0rgensen method to these clouds using our modified S/N requirements; namely, we used a S/N 
> 5 for at all IRAC bands in place of using a S/N > 5 at 70 /xm. 

We find the results from our technique are in reasonable agreement with those obtained from 
the J0rgensen and Enoch methods on the same data. There are some discrepancies and any 
given core could be classified differently between the various methods. For example, in Perseus, 
we identified 46 cores as protostellar, but only ~ 84% of these cores are similarly identified as 
protostellar with either the J0rgensen method or the Enoch method. Since the color selection 
criteria and coincidence criteria between these three methods are different, some discrepancies can 
be expected. 

In using intrinsic core FWHMs, the Enoch method is more prone to identify cores in clustered 
regions as protostellar, since there are generally more YSOs in these regions and the core sizes can 
be quite large (thereby, considering a greater coincidence area). For both Ophiuchus and Perseus, 
the Enoch method classified the most objects as protostellar. In using a fixed angular scale (ie., 
15") the J0rgensen method will cover a very different size scale at 125 pc (distance to Ophiuchus) 
compared to 450 pc (distance to Orion). As well, the SLC defined cores to have a minimum radius 
of 9.6" (based on a minimum of eight 6" pixels), therefore, depending on the intrinsic size of a 
core population and the distance (and the degree in which cores are unresolved and blended), the 
observed size may be smaller than a 15" search radius. For Taurus and Orion, the J0rgensen 
method identified many more protostars than our method, which could be related to the projected 
core size. Indeed, the average starless core size for Taurus and Orion (28.5" and 32.5", respectively) 
are smaller than the average size for Ophiuchus and Perseus (37.3" and 34.1", respectively), albeit, 
these differences are small. For Serpens, there are not enough cores to compare well the three 
methods. 

In comparison, we used a size based on a fractional flux to determine our coincidence criteria. 
This selection method considers both core shape and size, such that a smaller intrinsic core pop- 
ulation or more distant cloud will not be biased. In clustered regions, however, the submillimeter 
flux contours become complicated, so we are not fully capable of identifying these cores. As well, 
these clustered environments may have several bright sources that obscure or confuse the infrared 
emission towards nearby submillimeter cores (e.g., OMC-1, Peterson & Megeath 2008). There- 



^AU infrared sources in the c2d catalogue are labeled according to their observed SED. Designations oft;en asso- 
ciated with protostars include "YSOc" and "red". See the Final Delivery Document, Evans et al. (2007), for more 
details. 
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fore, identifying protostellar cores in clustered environments is more difficult, and we caution that 
extreme cases, like OMC-1, might require special consideration (see §4.1). 

To eliminate the problem of chance coincidences (in both clustered and isolated environments) , 
we applied a series of color criteria to eliminate obvious contaminants. Since we are interested in 
embedded sources, we select objects with strong MIPS detections (CCl and CC2) and colors that 
resemble YSOs (CC4a or CC4b). In addition, we exclude galaxy candidates based on colors (CCS) 
rather than magnitude, since the latter will vary with cloud extinction, thereby possibly removing 
deeply embedded protostars. Our criteria made no initial assumptions regarding c2d designation, 
and yet most of our protostellar sources in Ophiuchus, Perseus, and Serpens have c2d designations 
often associated with YSOs, suggesting that our method is capable of identifying protostellar sources 
without initially relying on c2d designations. This consistency allows us to apply our technique to 
data observed differently from the c2d clouds (Taurus and Orion) that have no such designations. 
There are, however, two protostellar cores with c2d designations of "galaxy candidates" in Perseus 
and one each in Ophiuchus and Serpens. 

It is important to note that the J0rgensen method and the Enoch method were developed 
based on different data sets. The J0rgensen method employed non-SLC SCUBA maps with a pixel 
resolution of 3" and a beam angular resolution ~ 15", whereas the SLC has a pixel resolution of 6" 
and a smoothed beam angular resolution ~ 23" at 850 ^m. Thus, the 15" distance limit used in the 
J0rgensen method may be inaccurate for the SLC data. The Enoch method was developed using 
cores identified at 1.1 mm (not 850 /nm) with Bolocam (40" FWHM resolution) for their respective 
analyses, and so the core numbers, locations and sizes from each sample will differ. 

3.5. Comparison to Known Class Sources 

Proebrich (2005) compiled a database of young protostars from the literature, and used all 
available photometry between 1 fim and 3.5 mm to build SEDs and re-classify the sources under 
a consistent system. This sample includes two objects in Ophiuchus and Taurus, ten objects in 
Perseus, and eighteen in Orion. For the most part, we identified objects in common as protostellar. 
Some protostars listed by Proebrich (2005), however, have no corresponding submillimeter cores in 
the Fundamental SLC. For example, the known protostar VLA 1623, which is located ~ 0.5' south- 
west of the bright and very crowded Oph A filament in Ophiuchus, does not have a submillimeter 
core in the Fundamental SLC. This core is listed in the Extended SLC, which we do not consider 
(see §2.3.1). 

In addition, some known young protostars were misidentified as starless cores due to no or 
poor infrared detections with Spitzer. For example, two protostars in Perseus, SVS 13B and NGC 
1333 IRAS 4A, had negative bandfilled IRAC fluxes and thus, failed our color criteria. For the 
most part, there is good agreement between our classification and the compilation from Proebrich 
(2005), thus our technique is capable of identifying protostellar cores, provided that Spitzer was 



-16- 



able to detect the source. Indeed, we are limited in our ability to classify our SCUBA cores by the 
areal coverage of Spitzer in relation to the SCUBA observations. To be as unbiased as possible, 
we did not alter our lists to include known protostellar cores that were labeled starless with our 
technique. We will consider separately the samples of starless and protostellar cores in the next 
section. 



4. Discussion 

One of our main goals is to examine CMFs for a variety of nearby clouds. In the previous 
section, we applied our own classification technique to identify which cores were starless and which 
were protostellar. To produce the CMFs, we must first estimate the mass contained within the 
cores. The mass of a SCUBA object can be estimated from the Planck function, 850 /im flux, 
opacity, and distance. For example, using the formula (Johnstone et al. 2000b): 



Mclump — 0.19 5'850 
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where S'gso is measured in Jy. For all clouds, we assumed a constant dust opacity, Kgso = 0.01 
cm^ g~^, though K850 can vary by a factor of two (Henning et al. 1995). Our cloud distances and 
dust temperatures are listed in Tables 1 and 2, respectively. We also assumed a single uniform 
temperature for all cores in a given cloud. Temperature distributions across cloud cores are not yet 
well known. Note that dust temperatures in starless cores should decrease from core edge to core 
centre as external heating from the interstellar radiation field is increasingly damped (Di Francesco 
et al. 2007). Protostellar cores, however, may have internal heating which will affect the radial 
temperature variations, such as a warm centre, cool middle and warm edges. 

Tables 10 and 11 in Appendix A list the starless and protostellar cores in Perseus identified 
using our classification technique. Each table lists core names (positions), peak fluxes, total fluxes, 
masses, and effective radii. The quoted total fluxes and effective radii are the "alternative flux" and 
"alternative radius", respectively (see §2.3.1). For complete tables of all 5 clouds, see the electronic 
edition. 



4.1. Starless CMFs 

Figure 4 shows the observed starless CMFs for the Ophiuchus, Taurus, Perseus, Serpens, and 
Orion clouds. These starless cores were identified using our classification scheme (see §3.2 and 
§3.3). To estimate errors in each mass bin, we varied the temperature in steps of ±1 K up to a 
30% uncertainty in T^, and found the standard deviation across the mass bins generated from the 
CMFs at these temperatures. Wc found such AT standard deviations were similar in magnitude to 
errors expected from Poisson statistics. In each panel of Figure 4, a power-law relationship with a 
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Salpeter slope (—1.35, see §4.1.1) is also shown as a dotted line and this seems to trace reasonably 
well the higher mass ends of the CMFs. 

Since we use the same classification scheme for all five clouds, differences in the CMFs should 
reflect differences between core populations in the clouds. For example, Orion contains more massive 
starless cores than Ophiuchus or Taurus, though many of the highest mass cores in Orion could be 
misidentifications associated with OMC-1 (see below). Conversely, Ophiuchus and Taurus contain 
some very low-mass cores. This could be related to cloud distance (see Table 1), as Ophiuchus (125 
pc) and Taurus (140 pc) are much closer than the other clouds, but a similar difference is not seen 
between Perseus (250 pc) and Orion (450 pc). 

In general (but excluding Serpens), the CMFs have a similar "lognormal" appearance, but 
some differences between them remain. Perseus shows the narrowest distribution. Orion has the 
widest distribution, due to the extreme high-mass extent. These high mass "cores" have masses on 
order of ~ lO^^"' Mq and may be larger-scale structures than cores. Recall that wc have defined 
a core as a dense, compact structure that could form a single star or a small stellar system. For 
Ophiuchus, Taurus, Perseus, and Orion, the CMF peaks are at roughly 0.1 Mq, 0.2 Mq, 1.6 Mq, 
and 1.0 M0, respectively. The Serpens CMF has a very small sample, and we cannot make the 
same conclusions or comparisons. 

As discussed before in §3.4, the clustered region, OMC-1, in Orion is a very complex region. 
OMC-1 contains bright infrared nebulosity which lowers the sensitivity to protostars and could 
result in cores misidcntificd as starless. In addition, OMC-1 is exposed to intense UV radiation due 
to nearby OB associations (Peterson & Megeath 2008). Therefore, it is also possible that we are 
under-estimating the temperature (20 K) in OMC-1. Indeed, increasing the temperature of these 
cores by a factor of 2 reduces the core masses by a factor of 2.5. To understand the impact the OMC- 
1 cores have on the Orion CMF, Figure 5 shows the Orion starless CMF after the cores specifically 
associated with OMC-1 (12 starless cores) are removed. Excluding these objects effectively removed 
the highest mass cores (M > 15 Mq) and made the CMF thinner, suggesting that the cores in 
OMC-1 may be different from the rest of Orion. To remain unbiased to our classification technique, 
however, we will mainly consider the total Orion CMF for the following analysis. We will also briefly 
compare the differences with the OMC-1 cores removed for the discussion of starless CMFs only. 

It is important to note that SCUBA (with a limited detection sensitivity) and Clumpfind 
(using a threshold of 90 mJy beam^^) are more sensitive to concentrated, bright cores. As such, 
our CMFs are incomplete for cores of low surface brightness. More sensitive instruments (i.e., 
SCUBA-2) will locate better these low surface brightness objects in the future, but we are unable 
to include these objects at present. For the most part, low surface brightness cores are expected to 
occupy the low-mass regime of the CMF, though there can be cases where high-mass "cores" are 
too diffuse to detect. The high- mass ends of our CMFs have been comparatively well sampled and 
these data can be compared best to the IMF. 
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4.1.1. Relation to the IMF 

Using an observed luminosity function for nearby stars, Salpeter (1955) found that the stehar 
mass distribution (or IMF) seemed to obey a power-law relation, ^(m) oc m~^-^^ for 0.4 Mq < M 
< 10 M0. Here, ^(m) is defined as dN = ^(m)dlogm, where dN is the number of stars of mass m 
lying between m and m + dm (Warner 1961). Thus, 



dN 
dlogm 

dN 
dm 
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In Figure 4, the Salpeter power-law slope appears to agree with the cloud CMFs at higher masses. 
To test this relationship more quantitatively, we calculated the ordinary least squares slope to the 
high-mass end of each CMF weighting each bin by its uncertainties. Since quantifying a high-mass 
end is uncertain, we determined several best fit slopes using different mass ranges with at least 4 
mass bins to calculate the best fit slope. Figure 6 shows the resulting best fits (solid lines) for one 
choice of mass range (dotted lines). For Orion, there is a slight increase in cores at very high masses 
(~ 10^~^ M0) mainly associated with OMC-1 (see Figure 5) and since these objects may not be 
cores, we do not consider these massive objects in our best fit slopes in Orion. Wc also calculated 
the best fit slopes for the Orion CMF with the cores towards OMC-1 removed using similar mass 
ranges. Table 7 gives the average ordinary least squares best fit for the well-sampled clouds and 
the best fit for Orion without the 12 starless cores from OMC-1. The Serpens CMF had too few 
mass bins to calculate a best fit slope. 

The slopes in Table 7 hint at some differences, though all agree with the Salpeter power-law 
slope, within ~ 2a errors. Indeed, some clouds have CMFs that are more Salpeter-like than other 
clouds. The errors were derived from the best fit and reflect the shape of the CMFs. For example, 
a CMF that appears curved (e.g., see Orion) will have a larger uncertainty associated with a linear 
slope than a CMF that follows a more linear drop-off (e.g., see Taurus). For the Orion CMF without 
the OMC-1 cores (see Figure 5), the shape is less curved and the errors more similar to what is 
found for Ophiuchus and Perseus (see also, Table 7). The slope obtained for the complete Orion 
CMF (-1.69±0.72) flattens considerably when the cores from OMC-1 are removed (-0.93±0.18). 

In addition to producing ordinary least squares fits, we also used the Kolmogorov-Smirnov (KS) 
Test on the high-mass end of our CMFs. We generated random core masses following a Salpeter 
distribution within the mass ranges listed in Table 7 and determined the likelihood that our observed 
CMFs were drawn from the same distribution as the random Salpeter sample. Unfortunately, the 
likelihood varied significantly with mass range, suggesting that our samples are too small for the 
KS test. A larger, more sensitive set of observations would significantly improve the comparison 
with the Salpeter IMF. 
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4.1.2. Trends with the CMFs 

Our five clouds represent a variety of environments and have very different properties (see §2). 
Any trends between these properties with the starless CMFs could reveal information about star 
formation across these environments. We compared the CMF peak mass and the best fit slope 

against cloud distance, cloud mass, cloud extinction peak, core line widths, core temperature, core 
extinction peak, and number of cores in the sample. We found very few correlations, in spite of our 
range of cloud properties. Figure 7 shows four of the fourteen trends that we examined. 

Figure 7a shows no correlation between core temperature and peak mass, whereas Figure 
7b (core temperature with best fit slope) shows the strongest case for a potential trend. The 
correlation in Figure 7b, however, depends greatly on the slope of Orion, which has the largest 
measured uncertainty. If you consider the slope obtained for Orion excluding the cores towards 
OMC-1 (—0.93 lb 0.18), this trend disappears. A trend between core temperature and best fit 
slope would imply that there is no universal power-law distribution for CMFs, unless starless cores 
exist at the same temperature. A more extensive study with more clouds is necessary to test this 
relationship (note, that our core temperatures are not particularly well-defined - see §2.2). Figure 
7c shows a possible trend between cloud distance and CMF peak mass. A slight positive relation 
between these properties may reflect our imposed threshold of 90 mJy beam~^. From Equation 3, 
mass should vary with distance and temperature, however there is no strong correlation between 
peak mass and temperature (Figure 7a). The weak trend in Figure 7c could also result from 
resolution problems, such as numerous lower mass cores appearing as a single higher mass cores in 
low resolution distant clouds. This would suggest, however, a shallower best fit slope with distance, 
which we do not observe in our furthest cloud, Orion, unless we exclude the OMC-1 cores. Figure 
7d shows no correlation between best fit slope and the number of cores. 



4.1.3. CMFs Using Different Classification Methods 

We applied two additional techniques in classifying a core as protostellar or starless: the 
J0rgensen method and the Enoch method. The J0rgensen method is outlined in J06, J07, and JOS, 

and used Spitzer c2d data and non-SLC SCUBA observations for Perseus (J06) and Ophiuchus 
(J07), whereas the Enoch method, outlined in EOS and E09, used Spitzer c2d data and Bolocam 
1.1 mm observations for Ophiuchus, Perseus, and Serpens. 

Figures S and 9 compare our starless CMFs for Ophiuchus and Perseus, respectively, with the 
starless CMFs obtained from classifying cores with the J0rgensen and Enoch methods. Although 
there are relatively few cores in general with large masses, the high-mass end of the CMF appears 
to be vary slightly with the classification technique. Table S compares the number of starless high 
mass cores in Ophiuchus and Perseus using the three techniques. We considered one mass range 
from our best fit results (see Table 7) to define the high-mass regime in each case. The Enoch 
method tends to classify the most massive cores as protostellar, resulting in the fewest that are 
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identified as starless. Indeed, a similar result by EOS led them to conclude that higher mass cores 
evolve more quickly than lower mass cores. All three methods, however, still produce CMFs that 
agree with the Salpetcr power-law distribution (dotted lines in Figures 8 and 9). The low-mass 
ends also show some variations (especially in Ophiuchus), but due to incompleteness in this regime, 
we cannot determine the significance of these differences. 

4.2. Protostellar CMFs 

In identifying starless cores with our classification technique, we also identify the protostellar 
core population. Figure 10 shows the protostellar CMFs for each cloud and their relation to 
a Salpctcr-Iikc slope. The protostellar cores in Figure 10 were identified using oTir classification 
technique. We initially assumed the same temperatures as with the starless CMFs and uncertainties 
were measured in the same manner as the starless CMFs. 

The protostellar CMFs do not show quite the same "lognormal" shape as the starless CMFs, 
though there are fewer numbers of protostellar cores than starless cores. Also, all protostellar 
CMFs seem to be systematically shifted to slightly higher masses than the starless CMFs. The 
Ophiuchus and Taurus protostellar CMFs peak at higher masses (~ 0.5 Mq for both) than their 
starless distributions (~ 0.2 Mq for both). The Perseus protostellar CMF peaks at ~ 2.5 Mq 
and the Orion protostellar CMF peaks at ~ 6.3 Mq, which are both higher than the peak for the 
starless CMFs (-- 1.0 Mq). 

Figure 10 seems to suggest that protostellar cores are intrinsically more massive than starless 
cores, but we assumed they have identical dust temperatures to the starless cores. Protostellar 
cores, however, may be heated internally and thus, have slightly higher temperatures (e.g., see 
EOS). A higher protostellar core temperature will reduce the core mass and possibly remove the 
apparent discrepancy between the starless and protostellar CMFs. 

We tested the effect of higher temperatures on the protostellar CMFs by raising the protostellar 
core temperatures in steps of AT^ = 1 K. We determined our final protostellar core temperature 
based on fitting the width and peak of the protostellar CMF to the starless CMF of a given cloud, 
assuming that these widths and peaks are similar and any envelope mass loss around a protostar is 
insignificant. In the end, we found good fit temperatures for the protostellar cores to be 20 it 2 K 
for Ophiuchus, 17 it 2 for Taurus, 15 it 1 K for Perseus, and 59 it 3 K for the entire Orion CMFs 
(Serpens has too few cores to properly match the width and peak of the CMF distributions). Thus, 
the protostellar core temperatures in Ophiuchus, Taurus, and Perseus should increase by a factor 
of ~ 1.3 and the temperatures in Orion should increase by a factor of ~ 3. Figure 11 illustrates the 
protostellar CMFs at their new temperatures in relation to the starless CMFs. The protostellar 
CMFs in Figure 11 show no overabundance in protostellar cores at higher masses relative to starless 
cores at the new temperatures, suggesting that higher mass cores do not necessarily evolve more 
quickly than lower mass cores, assuming insignificant mass loss. 
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The Orion protostellar CMF has a substantial temperature increase (59 K versus 20 K) com- 
pared to the other clouds, which may be related to "hot" cores often associated with massive star 
forming regions (Kurtz et al. 2000). Hot cores are very dense and can have temperatures of > 50 
K (Cesaroni et al. 1994). Figure 5 shows the Orion starless CMF with the high-mass star forming 
OMC-1 region removed. The resulting starless CMF is considerably more "lognormal" and consis- 
tent with the other clouds. This suggests that the cores in OMC-1 may be of unusually high-mass 
or at a higher temperature. 

For the rest of our analysis, we will compare starless and protostellar core properties using the 
same dust temperatures (see Table 2) and not the increased temperatures shown in Figure 11. In 
addition, for the rest of our analysis, we will consider the Orion core population with OMC-1 as 
opposed to removing the 12 starless cores in OMC-1 from our sample. 

4.3. Core Environments 

Figure 12 compares the masses and the observed radii of the cores in our five clouds. The 
observed radii are the true radii convolved with the SLC beam and then truncated according to the 
90 mJy beam~^ contour threshold. The SLC beam is approximated as a Gaussian with a FWHM 
of 22.9" (Di Francesco et al. 2008). Due to this finite resolution, all cores will have smaller true sizes 
than observed, with less extended cores more affected than larger ones. We have not deconvolved 
the sizes of our cores in Figure 12. 

For all the clouds in our sample, core mass and radius are well correlated. Figure 12 also 
includes the sensitivity limit (dot-dashed line) from the 90 mJy beam~^ threshold. Although the 
scatter is quite restricted in one direction due to sensitivity, we note that at higher masses, core 
mass begins to rise faster than size in Figure 12. Indeed, Figure 12 illustrates our completeness 
level. At the lowest masses, we are very close to the limit of our ability to detect cores (closer to 
our sensitivity limit). For Serpens, the protostellar cores appear to have higher masses and larger 
sizes than the starless cores. This trend, however, may be a reflection of the small sample, since 
the other clouds do not show the same results. 

Figure 13 gives the Ay distributions for the five clouds, separated by protostellar (solid lines) 
and starless (dashed lines) cores. The histograms are all binned to A Ay = 4 to ensure each bin is 
well populated. Within a given cloud, the distributions arc broadly similar. For example, both the 
protostellar and starless core Ay distributions for Ophiuchus peak at Ay ^ 25 magnitudes. This 
peak is at significantly higher Ay than the peaks in the other clouds, though this could be related 
to distance. For more distant clouds, we are likely less sensitive to higher extinctions due to lower 
linear resolution. Serpens is the only exception to the similar peaks in Ay between protostellar 
and starless core populations, though that is likely due to its small sample size. 

In Figures 14 and 15 we compare extinction with core mass and core size, respectively. For 
Ophiuchus, Taurus, Perseus, and Orion, we found no strong trends between these distributions. 
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Serpens may show trends of increasing mass and size with extinction for starless to protostellar 
cores, though the other clouds with larger samples do not reveal the same trends. For the other 
clouds, the protostellar and starless cores have wide ranges of mass and size over a wide range of 
extinctions. This finding is counter to the expectation that cores at higher extinction could be 
at higher pressures and hence may be smaller in size (Johnstone et al. 2004). While there is a 
correlation between core size and mass (see Figure 12) for all five clouds, no obvious relations of 
these quantities with extinction are seen given a good sample size. 

4.4. Predicted CMFs 

SCUBA was limited in its ability to sample the dense core populations in nearby star-forming 
clouds and typically observed known regions of active star formation (e.g., L1688, Johnstone et al. 
2000b), so very large fractions of these clouds, particularly at low extinctions, remain unmapped 
(see Figure 1). As a result, the CMFs derived from SLC map data (see §4.1) do not represent all 
cores within their respective clouds. 

SCUBA-2 will be the successor continuum mapping instrument to SCUBA on the JCMT. The 
SCUBA-2 Gould Belt Legacy Survey (GBLS) will map 15 star-forming molecular clouds within 
500 pc with excellent sensitivity to detect core down to substellar masses (Ward-Thompson et al. 
2007b). These data will allow us to obtain robust CMFs, construct a less ambiguous approach to 

classifying cores, constrain the lifetimes associated with star formation, and solidify any possible 
similarities between CMFs and the IMF. With superior sensitivity and efficiency, SCUBA-2 will be 
able to map nearby clouds with improved flux sensitivity and speed, allowing for more complete 
samples of their core populations. 

We wish to extrapolate our current CMFs to predict the CMFs we will obtain using SCUBA-2. 
Predicted CMFs can be created by assuming the observed CMFs trace the incidences of cores for 
the unobserved regions of each cloud. Of course, since SCUBA-2 will have improved sensitivity, 
lower surface brightness cores will be more easily detected. As well, our sample is limited to cores 
with minimum peak fluxes of 150 mJy beam~^ and size limits of 90 mJy beam~^ (following our 
original criteria for identifying cores; see §2.4). Thus, we caution that our predicted CMFs are only 
applicable in the higher mass regime where our observed data are more complete. 

The predicted starless CMFs will be extrapolations of the observed starless CMFs over the 
unobserved regions of each molecular cloud. We considered two techniques to create our predicted 
CMFs. First, an ^dy-independent extrapolation where we gave each core equal weight and extrap- 
olated the observed starless CMFs using the ratio of total cloud area (for Ay > 3) to cloud area 
observed by SCUBA. Second, an ^y-dependent extrapolation, where we extrapolated core inci- 
dences over small extinction ranges using the ratio of total cloud area over ranges in Ay to cloud 
area observed by SCUBA over the same ranges of Ay- Equation 6 outlines the Ay -dependent 
extrapolation. We chose AAy = 4 to ensure that the extinction ranges are well populated. For 
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a given Ay range, the predicted number of objects {NAy,pred) is given by the observed number in 
that range (-A^i^ q^s) multiphed by the ratio of entire cloud area to the observed cloud area at that 
Ay range. 

.cloud /n\ 

J-^Av,pred = J^^Av,obs7^ (.Oj 

Ay, scuba 

where ^ Ay, cloud is the area of the whole given cloud that falls within a given Ay range and ^ Ay, scuba 
is the area within the same extinction range observed by SCUBA. 

The fraction of each cloud observed by SCUBA was different from cloud to cloud. Thus, taking 
into account the area of observations can help us compare the results between clouds. Each SLC 
map is 1.2' x 1.2' in extent and several overlap with other maps (Di Francesco et al. 2008). Table 
3 lists the total areas observed by SCUBA for each cloud. 

One additional complication is that the Scorpius and Monoceros clouds have projected locations 
near Ophiuchus and Orion, respectively (see Figure 1). As such, the 2MASS extinction maps of 
Ophiuchus and Orion contain sections of Upper Scorpius and Monoceros R2. For Serpens, we 
consider the entire 2MASS map since the cloud boundary is not clear as with the other clouds. 
Since our core counts with Serpens are quite low, this will not greatly affect our overall results. 



4.4- 1- Predictions 

The predicted CMFs are shown in Figure 16 for each cloud, comparing the observed starless 
CMF (dotted histogram) to the Ay-dependent extrapolated (solid histogram) and Av^-independent 

extrapolated (dashed histogram) starless CMFs. Overall, the predicted starless CMFs are generally 
similar in shape, even in the ^ly-dependent case, but have large differences in number. 

Table 9 lists the numbers of observed starless cores and the numbers of predicted starless cores 
using the above assumptions. The quoted errors in our predicted core counts were determined from a 
quadrature sum of the Vn uncertainties in each bin of the observed CMF scaled upwards according 
to the ratio of unobserved to observed area. The errors given in Table 9 show that the difference 
between the predicted numbers and the observed numbers are significant. Our analysis, however, 
is biased against low-mass objects, particularly diffuse sources, which were more difficult to detect 
or keep within our core definition limits. Since we cannot accurately probe the low-mass end of the 
observed CMFs, our predicted CMFs are also incomplete in this regime and the numbers quoted in 
Table 9 are lower limits to the actual expected core statistics using these methods. Nevertheless, 
we predict increases in the total numbers of starless cores by factors of at least 6 - 11 in future 
surveys. 

There is a significant increase in the number of predicted cores obtained from using the 
Ay-independent extrapolation over the numbers from the Ay-dependent case. Since the Ay- 
independent method predicts so many more cores over the Ay-dependent method and the majority 
of the cloud exists at lower extinction levels, this increase suggests that cores observed by SCUBA 
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tended to be located in regions of relatively high extinction. Indeed, SCUBA mainly focused on 
the very high extinction regions of each cloud (see Figure 1). 

4-4- ^- Gaussian Replacements 

The predicted CMFs in Figure 16 were made assuming we accurately knew the mass of each core 
based on flux and temperature. This assumption, however, may be unwarranted since uncertainties 
from the instruments, temperature gradients, distance measurements, and cloud opacities will all 
affect the estimated core mass. We attempted to simulate this uncertainty by replacing core 
populations in every mass bin of the observed CMFs with a Gaussian of the same area. We set the 
FWHMs of our Gaussian replacements to be the width of our initial mass bins {log M/Mq = 0.2) 
and with these replacements, we computed predicted CMFs with the same assumptions as before. 

Figure 17 shows the predicted starless CMFs assuming an Ay-dependent extrapolation and 
Figure 18 shows the predicted starless CMFs for the Ay-independent extrapolation. Included in 
both figures are the cumulative Gaussian replacements for each of the mass bins and a Salpeter 
power-law slope for comparison. The Gaussian replacements show a low-mass "tail" due to pro- 
jecting the Gaussians from a linear plane to a logarithmic plane. The errors in the bin populations 
were determined by varying the temperature with the observed starless CMFs as before. The slopes 
of the predicted histogram distributions and the predicted Gaussian replacement distributions all 
agree with the Salpeter power-law, within errors. 

5. Conclusions 

We have used data from SCUBA, Spitzer, and 2MASS to produce relatively unbiased starless 
and protostcllar CMFs using a new classification technique. We applied this technique to the 
Ophiuchus, Taurus, Perseus, Serpens, and Orion molecular clouds. To our knowledge, this work is 
the most extensive study of the starless core populations in nearby molecular clouds to date, since 
we have examined 5 unique clouds in a uniform manner where previous studies have focused on 
smaller regions or fewer clouds (e.g., J07; Nutter &; Ward-Thompson 2007; EOS; Nutter et al. 2008) 
with different core identification schemes. 

It is important to consider that our core populations are incomplete. We were unable to 
detect well cores of low surface brightness (diffuse) with SCUBA. These cores will be more easily 
observed with more sensitive instruments, such as SCUBA-2. We expect that the low-mass end 
of our CMFs will be most affected by the improved sensitivity of SCUBA-2. The high-mass end, 
however, should be relatively complete. As well, we have developed a self-consistent technique 
to analyze these five clouds and determine which cores were starless and which were protostellar. 
Our technique identifies YSO-like infrared sources through strong infrared emission and red colors 
associated with embedded protostars. We then used the 850 fim flux at the position of each of 
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these YSO-like sources to determine which were associated with our submilhmeter cores. With our 
classification technique, we developed uniform samples of cores for each cloud. 

Assuming a constant dust temperature for all cores in a given cloud, our starless CMF best fit 
slopes were -1.26 ±0.20, -1.22 ±0.06, -0.95 ±0.20, and -1.67 ±0.72 for the Ophiuchus, Taurus, 
Perseus, and Orion starless CMFs, respectively. Removing the OMC-1 cores from Orion flattened 
the slope to —0.93 ± 0.18. These slopes all agree with the Salpeter power-law (Salpeter = —1.35) 
within errors, though it is possible that the OMC-1 region represents a population of cores with 
very different properties with respect to the rest of the Orion cloud. We were unable to calculate 
a best fit slope for Serpens due to the small sample size. 

Our starless CMFs showed similar log-normal shapes, though the range of masses and peak 
in the distribution varied with each cloud. The starless CMF for Orion also included an unusual 
population increase at very high masses due to cores in OMC-1, which could be attributed to its 
cores having a more varied temperature or to difficulties in distinguishing between starless and 
protostellar cores with the infrared data. Nevertheless, with our self-consistent technique, our four 
well sampled clouds all agree with a Salpeter power-law distribution within 2cr errors (except Orion 
without OMC-1 which is ~ 2.3cr). To better understand the shape of these CMFs and their relation 
to the IMF, we need larger samples to minimize the uncertainties. 

We also tested our classification scheme against two others (the J0rgensen and Enoch methods) 
and found that the fraction of high-mass cores identified as starless or protostellar varied slightly 
with the technique used for classification. In addition, we also examined trends in the starless CMF 
peak masses or best fit slopes with cloud properties, and we found a potential relationship between 
best fit slope and core temperature, suggesting that a single power-law relationship for CMFs may 
not be ideal. This correlation, however, depends greatly on the slope of the Orion CMF, which has 
the largest uncertainty. Indeed, this trend disappears if we consider the Orion CMF with the cores 
towards OMC-1 removed. Since the CMFs all agree with a Salpeter distribution within errors, 
more data at better resolution and good sensitivity are necessary to understand these relations. 

Assuming the same dust temperature as the starless cores, the protostellar CMFs were sys- 
tematically shifted to slightly higher masses. We increased the protostellar core temperatures to fit 
the width and peak of the protostellar CMFs with their starless CMFs and found the protostellar 
core temperatures should increase by a factor of ~ 1.3 in Ophiuchus, Taurus, and Perseus, and a 
factor of ~ 3 in Orion over our original dust temperatures. Though these increases are plausible, 
radiative transfer modeling of protostellar cores is necessary to determine if protostellar heating 
has such a significant impact on the average core temperature. 

Finally, we used simple assumptions to predict the core populations for the unobserved regions 
of each cloud. We developed two sets of predicted CMFs, an ^y-dependent extrapolation and 
an ^y-independent extrapolation, both which agreed with a Salpeter power-law distribution. We 
found that the ^y-independent extrapolation predicted a larger number of dense cores in each cloud 
than the ^y-dependent extrapolation had predicted. We hope to compare our predicted number 
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of cores to the results from surveys of entire cloud regions in the future (ie., with SCUBA-2). 
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APPENDIX 

A. Starless and Protostellar Core Lists 

Using the SLC and our 2MASS data, wc identified cores and their locations within the bulk 
cloud. Tables 10 and 11 give examples of the SLC and extinction information for cores in the 
Perseus cloud. For both tables, the core name reflects the J2000 position of the submillimeter peak 
and the mass is calculated using Equation 3 with the temperatures given in Table 2. The complete 
tables are available in the electronic edition. 
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Fig. 1. — 2MASS extinction maps with Funamental SLC data (contours). For Ophiuchus and 
Orion, lines are also shown to mark the rough boundary between the Scorpius and Monoceros 
clouds, respectively. For the color figure, see the electronic edition. 




[B.0]-[24] 



Fig. 2. — Color-color diagrams for the Spitzer objects in our sample. The diamonds show the 
infrared objects that met all our color criteria. The grey circles are objects that passed the high 
quality conditions (CCl and CC2) but with an associated 850 /im flux of 5*850 < 90 mJy beam~^. 
The crosses are objects with an associated 850 //m flux of S'sso > 90 mJy beam~^ but also met the 
colors of a star- forming PAH galaxy (CCS). The large concentration of grey points in Taurus at 
[8.0] — [24] ~ are due to foreground and background stars (the c2d catalogue used a pipeline that 
identified and removed many stars from the catalogue). As well, a number of infrared sources in 
Taurus were saturated, resulting in a "rollover" at [3.6] — [5.8] ~ 3. For the color figure, see the 
electronic edition. 
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Fig. 3. — Schematic of our classification technique for protosteUar cores. A subminimeter core is 
considered to be protostellar if a YSO candidate identified in §3.2 falls within the shaded region 
illustrated in this figure. The black contour represents the sensitivity boundary of the core (90 mJy 
beam~^) and the pale grey circle represents the peak intensity. The dark grey contour represents 
the maximum distance from the central peak for the protostar classification (as given by Equation 
2). 
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Fig. 4. — Starless CMFs for Ophiuchus, Taurus, Perseus, Serpens, and Orion with Salpeter power- 
law slopes dN/d{logm) oc m~^'^^. Uncertainties were determined from varying the temperature 
(see text). Starless cores were identified using our classification technique outlined in §3.2 and §3.3. 




Fig. 5. — Same as Orion in Figure 4 but with cores towards OMC-1 removed. We used coordinates of 
[5^35™00^ 5''35"*32^] in right ascension and [-5°29'00", -5°17'00"] in dechnation as our boundary 
for OMC-1. 
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Fig. 6. — Best fit slopes from ordinary least squares regression. The dotted lines indicate one 
mass range used to calculate the slope (see Table 7). The Serpens CMF had too few mass bins to 
calculate the best fit slope. 



-39- 



40 ^ ' ' ' I 

- A Ophiuchus 

35 - tK Taurus 

: O Perseus 

~_ • Serpens 

30 - □ Orion 



15 - 
10 - 
5 _ 



peak mass (Mq) 



A Ophiuchus 
H Taurus 
O Perseus 
• Serf:ens 
□ Orion 



-1.5 -2.0 

best fit slope 



A Ophiuchus 
M Taurus 
O Perseus 
• Serpens 
□ Orion 



A Ophiuchus 

Taurus 
O Perseus 
• Serpens 
□ Orion 



peol< moss (ivlg) 



-1.5 -2.0 
best fit slope 



Fig. 7. — Examples of comparisons between cloud or core properties with CMF properties. Peak 
masses were defined as the most populated mass bins in the starless CMFs. The best fit slopes are 
given in Table 7. Serpens had too few points to determine a linear best fit slope. 
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Fig. 8. — Comparison of the resulting CMFs from applying the EOS, J07 and our own core classifica- 
tion techniques for Ophiuchus. The dotted line represents a Salpeter power-law slope. Uncertainties 
were determined by varying the temperature by ~ 30%. 
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Fig. 9. — Same as Figure 8 but for Perseus. 
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Fig. 10. — Protostellar CMFs for all five clouds with Salpeter power- law slopes dN/dlogm oc 
j.^-1.35 Uncertainties were determined from varying the temperature up to 30%. We assume the 
same temperature as the starless cores. 



-43- 



1.0 



0.1 



10.0 



1.0 



0.1 



Ophiuchus 



Protostellor 
Storless 



Tourus 



Protostellor 
Storless 




Perseus 




Protostellor 
Storless 



Orion 



0.01 0.1 



1 

M/Mg 



10 100 1000 



0.01 



1000 



Fig. 11. — Protostellar (dotted) and starless (solid) CMFs for all five clouds. The temperatures 
for the protostellar CMFs have been adjusted to mateh the width and peak of the starless CMF. 
The temperatures used are 20 K, 17 K, 15 K, and 59 K for Ophiuehus, Taurus, Perseus, and Orion 
protostellar CMFs, respeetively. The respective starless core temperatures are 15 K, 13 K, 11 K, 
and 20 K. Serpens had too few cores to make this comparison. 
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Fig. 12. — Comparison of radius and mass for cores in the five clouds in our sample. Protostellar 
cores were classified using our technique outlined in §3.2 and §3.3. We assume a constant dust 
temperature (see Table 2) for cores within each cloud regardless of whether the cores are protostars 
or starless. The threshold sensitivity of 90 mJy beam~^ is shown as a dot-dashed line. 



-45- 




Fig. 13. — Distribution of extinction for each core in our five clouds. Cores were classified into 
starless or protostellar using the technique we developed (see §3.2 and §3.3). 
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Fig. 14. — Core extinction versus mass for starless and protostellar cores. Cores were classified 
using the technique outlined in §3.2 and §3.3. The dashed line illustrates the Ay = 3 limit we 
imposed on all cores. 
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Fig. 15.— 



Same as Figure 14 but showing core extinction versus size. 
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Fig. 16. — Observed starless CMFs are shown as dotted histograms for all five clouds with the 

predicted Ay-dependent starless CMFs as solid histograms and predicted ^y-indcpcndcnt starless 
CMFs as dashed histograms. Extinction ranges for the ^y-dependent starless CMFs of all clouds 
were AAy = 4. For Serpens, we considered the entire extinction map region to produce the 
extrapolated predicted CMFs (see text). 
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Fig. 17. — Predicted starless CMFs for all five clouds using the Ay-dependent extrapolation 
method. The uncertainties were measured by finding the standard deviation from changes in the 
temperature. The dashed lines represent a Salpeter slope. Also included is the final distribution 
obtained from replacing the counts in the mass histograms at each Ay range by a Gaussian (dotted 
curve). See text for more details. 
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Fig. 18.— 



Same as Figure 17 but for the Ay-independent extrapolation. 
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Table 1: Cloud Properties 



Cloud 


D (pc) 


Reference^ 


M(M0) 


Area (deg^) 


Ophiuchus 


125 


Enoch et al. 2009 


1.0 X 10^ 


67 


Taurus 


140 


Goldsmith et al. 2008 


3.9 X 10^ 


250 


Perseus 


250 


Enoch et al. 2009 


3.3 X 10'' 


66 


Serpens'' 


260 


Enoch et al. 2009 


1.9 X 10" 


12.3 


Orion 


450 


Peterson & Megeath 2008 


2.6 X 10^ 


147 



"Reference for our assumed distances. 

''Since the boundary between the Serpens and Auriga molecular clouds is unclear, we re-measured the 

mass of Serpens using the cloud boundary from the Hcrschcl Gould Belt survey (sec http://starformation- 
herschel.iap.fr/gouldbelt/gouldbelt_stage2.pdf), and we obtained a mass of 1.0 x 10'' Mq over 6.3 deg^. 

Table 2: Assumed Properties 



Cloud Td (K) Reference^ 

Ophiuchus 15 Priesen et al. 2009 

Taurus 13 Andre et al. 2000 

Perseus 11 Rosolowsky ct al. 2008 

Serpens 17 Schnee et al. 2005 

Orion 20 Wilson et al. 1999 



"Reference for our assumed temperatures. 
Table 3: Area Observed by Each Survey 



Cloud 


SCUBA 
(deg2) 


Spitzer 
(deg2) 


2MASS^ 
(deg2) 


Ophiuchus 


2.4 


6.6'' 


102 


Taurus 


0.94 


44b 


327 


Perseus 


2.7 


3.9'' 


99.8 


Serpens 


0.05 


0.85'' 


12.3 


Orion 


1.4 


9.95^= 


214 



"Areas of the entire 2MASS maps. For Ophiuchus and Orion, the 2MASS maps were edited to remove the Scorpius 
and Monoceros clouds, respectively (see Figure 1). 
''Area with both MIPS and IRAC data. 

■^Area with complete 4-band IRAC coverage. MIPS observations nearly covered the same regions. 
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Table 4: Summary of Cuts to the SCUBA Core Candidate List 



Cloud 


Initial^ 


Av>S 


.^850 > 0^ 


Visual 


Speak > 0.15 Jy beam""^ 


Ophiuchus 


151 


150 


150 


134 


124 


Taurus 


172 


162 


147 


117 


87 


Perseus 


246 


232 


231 


170 


143 


Serpens 


19 


19 


19 


15 


15 


Orion 


448 


436 


431 


392 


375 



"Number of submillimetcr core candidates within our 2MASS extinction maps. 

''All core candidates from the SLC were initially identified using a Clumpfind threshold of 3a, where a is the noise 
level in each map (see §2.3.1). We used the alternative flux, however, which measured the flux for each core again, 
but with a Clumpfind threshold of 3 x 30 m,Jy beam^^, where 30 m.Jy beam~^ represents the average noise in all the 
SLC maps. In low-noise maps, the peak 850 /am flux can be < 90 mjy beam^^, and thus, an alternative flux cannot 
be obtained. 



Table 5: Remaining Objects After Each Cut 



Cloud 


001-2*^ 


S850 > 90 


CC3 


CC4a-b 


Ophiuchus 


3474 


146 


116 


73 


Taurus 


91771 


167 


160 


33 


Perseus 


3429 


159 


133 


99 


Serpens 


1443 


34 


32 


25 


Orion 


18908 


803 


768 


226 



"CCl and CC2 for Ophiuchus, Perseus and Serpens included S/N > 5 for either 24 or 70 fun and that neither were 
bandfiUed. Orion and Taurus were not bandfiUed, and we used S/N > 5 for the 24 fim or for all IRAC bands. 
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Tal)l(^ (i: Comparison of Proloslcllar and Starless Core Numbers 



C loud 


A Lei hod 


Tt i ill 

rroloslc'itar 


C i 1 

blartcss 




J0rgensen 


25 


99 


Ophiuchus 


Enoch 


33 


91 




this work 


27 


97 




J0rgensen 


24 


63 


Taurus 


Enoch 








this \\"ork 


18 


G9 




J0rgensen 


42 


101 


Perseus 


Enoch 


49 


94 




this work 


46 


97 




J0rgensen 


8 


7 


Serpens 


Enoch 


7 


8 




tliis work 


11 


i 




J0rgensen 


109 


266 


Orion 


Enoch 








this work 


83 


292 



Table 7: Mean Best Fit Slope 



Cloud Slopes^ Mass Ranges^ (M©) 

Ophiuchus -1.26 ± 0.20 0.3 < M < 5; 0.5 < M < 5; 0.8 < M < 5 

Taurus -1.22 ± 0.06 0.3 < M < 3; 0.5 < M < 3 

Perseus -0.95 ± 0.20 1.0 < M < 20; 1.2 < M < 20; 1.6 < M < 20 

Orion -1.69 ± 0.72 1.4 < M < 32; 2.0 < M < 32; 3.2 < M < 32 



Orion (no OMC-1) -0.93 ± 0.18 1.4 < M < 20; 2.0 < M < 20; 3.2 < M < 20 



"The slopes quoted for the clouds were calculated from ordinary linear regression with each point weighted by the AT 
standard deviation method (sec text). The slope uncertainties come from the linear fitting algorithm. 

''Our best fit slopes were calculated using sixlin developed by Isobe et al. (1990), and we constrained the fits to have 
at least 4 points of reference. Serpens is not included because there were not enough mass bins to calculate a least 
squares fit. 
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Table 8: Comparison of High-Mass Starless Cores 



Cloud 


Method 


Mass Range 


Starless 




J0rgensen 


M > 0.5 M0 


33 


Opliiuchus 


Enoch 


M > 0.5 M0 


25 




this work 


M > 0.5 Mq 


27 




J0rgensen 


M > 1.2 Mq 


88 


Perseus 


Enoch 


M > 1.2 Mq 


80 




this work 


M > 1.2 Mq 


85 



Table 9: Predicted and Observed Starless Core Numbers 



Cloud 


Observed 


^y-dependent 


Ay-independent 


Ophiuchus 


97 


109 ± 13 


455 ± 46 


Taurus 


69 


467 ± 87 


1239 ± 149 


Perseus 


97 


165 ± 18 


232 ± 24 


Serpens 


4 


91 ±50 


776 ± 388 


Orion 


292 


2317 ± 182 


3476 ± 203 


Total 


559 


3149 ± 350 


6178 ±810 
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Table 10: Examples of Starless Cores^ 



Cloud 


Core Name 


j\v 




'-'850 


iVi 


^eff 






( Tonnn\ 
(^JzUUUj 


(mag) 


(Jy beam~^) 


( T-,r"\ 

(Jy) 


(J^o) 


(PC) 






J162538.4- 242238 


10.52 


0.15 


0.03 


0.0073 


4.1 X 10' 


-3 


Ophiuchus 


J162542.8-241708 


11.77 


0.15 


0.30 


0.0733 


1.4 X 10" 


-2 




J162546.8-241714 


12.36 


0.15 


1.06 


0.2589 


2.5 X 10" 


-2 




J040443.3+261859 


8.16 


0.40 


0.60 


0.2355 


1.7 X 10- 


-2 


Taurus 


J040447.9+261923 


8.16 


0.34 


2.93 


1.1498 


3.6 X 10" 


-2 




J041828.6+282735 


16.91 


0.15 


0.02 


0.0079 


4.0 X 10- 


-3 




J032544.0+304026 


11.63 


0.34 


0.25 


0.4279 


2.0 X 10- 


-2 


Perseus 


J032703.2+301513 


5.10 


0.23 


0.83 


1.4207 


4.1 X 10" 


-2 




J032829.4+310956 


4.72 


0.21 


0.96 


1.6433 


4.5 X 10" 


-2 




J183001.4+010903 


14.66 


0.40 


0.70 


0.6035 


2.8 X 10" 


-2 


Serpens 


J183002.6+010827 


12.86 


0.42 


1.00 


0.8621 


3.5 X 10" 


-2 




J183005.0+011515 


14.16 


0.21 


0.24 


0.2069 


2.2 X 10- 


-2 




J051938.8-055144 


7.55 


0.15 


0.15 


0.3020 


3.4 X 10- 


-2 


Orion 


J051940.4-055150 


8.22 


0.17 


0.11 


0.2215 


2.9 X 10" 


-2 




J051948.1-055202 


9.71 


1.07 


3.39 


6.8254 


9.5 X 10" 


-2 



"This table is published in its entirety in the electronic edition. A portion is shown here for guidance regarding its 
form and content. All cores were classified using our technique. 
''Visual extinction measured towards the location of peak SCUBA 850 //m flux. 
°The peak SCUBA 850 fmi flux (at the position of the source). 

''Total SCUBA 850 /xm flux within the core, where the core was defined by Clumpfind with a threshold of 90 mjy 
beam~^. 

'^Mass derived from total SCUBA 850 irni flux using Equation 3 with the assumed temperatures listed in Table 2. 
■'^Tlic effective radius of the core determined by the Clumpfind area, A, where -Re// = -^^/tt for a Clumpfind threshold 
of 90 mJy beam~^. 
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Table 11: Examples of Protostellar Cores'^ 



Cloud 


Core Name 




<? ^ 


'-'850 


A /re 
iVl 










(mag) 


(Jy beam~^) 


(Jy) 




(PC) 






J162540. 5-243020 


14.10 


0.23 


1.72 


0.4201 


2.8 X 10" 


-2 


Ophiuchus 


J162551. 5-243256 


22.52 


0.15 


0.27 


0.0659 


1.2 X 10" 


-2 




J162610.4-242056 


29.68 


0.65 


2.40 


0.5861 


2.7 X 10" 


-2 




J041354.4+281129 


24.83 


0.19 


0.36 


0.1412 


1.6 X 10" 


-2 


Taurus 


J041426. 2+280601 


14.23 


1.03 


0.85 


0.3336 


1.4 X 10" 


-2 




J041828.6+282735 


15.07 


0.69 


1.52 


0.5965 


2.4 X 10" 


-3 




J032522.2+304514 


24.21 


1.41 


2.13 


3.6460 


3.9 X 10- 


-2 


Perseus 


J032536.1+304514 


10.11 


5.46 


10.09 


17.2712 


5.6 X 10" 


-2 




J041942.3+271337 


7.54 


2.23 


2.88 


4.9298 


3.7 X 10" 


-2 




J182948. 2+011639 


27.01 


4.16 


13.21 


11.3883 


6.2 X 10" 


-2 


Serpens 


J182949.8+011515 


25.54 


9.87 


15.17 


13.0708 


5.5 X 10" 


-2 




J182951.4+011633 


17.68 


1.76 


3.91 


3.3708 


3.8 X 10- 


-2 




J053443.8-054126 


6.85 


0.48 


0.51 


1.0268 


4.6 X 10- 


-2 


Orion 


J053453.0-054138 


7.07 


0.19 


0.49 


0.9866 


5.6 X 10" 


-2 




J053502.3-053756 


9.31 


1.01 


3.67 


7.3892 


7.8 X 10" 


-2 



"This table is published in its entirety in the electronic edition. A portion is shown here for guidance regarding its 
form and content. All cores were classified using our technique. 
''Visual extinction measured towards the location of peak SCUBA 850 //m flux. 
°The peak SCUBA 850 /um flux (at the position of the source). 

''Total SCUBA 850 /xm flux within the core, where the core was defined by Clumpfind with a threshold of 90 mjy 
beam~^. 

'^Mass derived from total SCUBA 850 /um flux using Equation 3 with the assumed temperatures listed in Table 2. 
■'^Tlic effective radius of the core determined by the Clumpfind area, A, where -Re// = -^^/tt for a Clumpfind threshold 
of 90 mJy beam~^. 



